Microscopic dissipation in a cohesionless granular jet impact 
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Sufficiently fine granular systems appear to exhibit continuum properties, though the precise 
continuum limit obtained can be vastly different depending on the particular system. We investigate 
the continuum limit of an unconfined, dense granular flow. To do this we use as a test system a 
two-dimensional dense cohesionless granular jet impinging upon a target. We simulate this via a 
timestep driven hard sphere method, and apply a mean-field theoretical approach to connect the 
macroscopic flow with the microscopic material parameters of the grains. We observe that the flow 
separates into a cone with an interior cone angle determined by the conservation of momentum and 
the dissipation of energy. From the cone angle we extract a dimensionless quantity A — B that 
characterizes the flow. We find that this quantity depends both on whether or not a deadzone — 
a stationary region near the target — is present, and on the value of the coefficient of dynamic 
friction. We present a theory for the scaling oi A — B with the coefficient of friction that suggests 
that dissipation is primarily a perturbative effect in this flow, rather than the source of qualitatively 
different behavior. 



PACS numbers: 45.70.Mg, 47.57.Gc, 02.70.Ns 

In surprising experiments by Cheng, et al[l], it has 
been shown that a dense cohesionless granular jet impact 
can behave like a water flow of the same geometry. The 
experiment consists of a jet — either granular or liquid 
— being projected at some speed at a cylindrical target. 
The fluid forms a cone upon leaving the target, and the 
interior angle of the cone is measured (Fig. [l]) . The gran- 
ular flow, despite being cohesionless, is observed to form 
a coUimated cone when leaving the target. The interior 
angle of this cone quantitatively agrees with the cone an- 
gles observed when a water jet impinges upon a target of 
the same geometry [2]. Furthermore, the Reynolds num- 
ber of the water jet is very large, suggesting that dissipa- 
tion is not involved in determining the angle of the cone 
in the water flow. The comparison between the water 
flow and granular flow is therefore even more surprising, 
as the granular flow has strong dissipation from inelas- 
ticity and friction. This suggests somehow that despite 
the presence of large dissipation in the granular flow, that 
dissipation is not strongly influencing the bulk properties 
of the flow. 

This is in contrast to the majority of granular sys- 
tems, in which the dissipation is critical in determining 
their behaviors. For example, granular shear flows ex- 
hibit ordering depending on their inelasticity [31 |5j, and 
continuum theory of granular flows has a singularity with 
respect to the dissipation's]. Even more surprising, flow 
past an obstacle - virtually the same geometry as the 
jet impact experiment - has very different phenomena 
when the granular flow is more dilute, that seem more 
in line with the dissipation-dominated picture of granu- 
lar flows. Granular flows impinging upon obstacles have 
been observed to form Shockwaves both in experiment f5^ 
and simulation 7 . The behavior of these shocks and the 
drag on immersed objects is seen to be dependent on in- 
elasticity, whereas the cone angle observed in the dense 



granular flow seems to be completely independent of the 
inelasticity of the material used in the experiments [T]. 

There is an additional aspect of the dense granular flow 
that makes its correspondence to water flows surprising: 
the granular jet forms a deadzone — a region in which 
particles are trapped in a static arrangement — near the 
target. This type of flow structure is seen in granular 
column collapse, in which a certain part of the column 
does not displace throughout its dvnamics J8l-fl0] . This 
is a manifestly granular behavior that does not appear 
anywhere in the corresponding liquid jet case. Changing 
the target shape influences the cone angle for water jTT], 
but adding this dead region in the granular flow does not 
seem to alter the agreement between the water flow cone 
angle (which lacks such an internal shape) , and the gran- 
ular flow cone angle. One might expect that the dead- 
zone would have a similar impact on the granular flow as 
a wedge would on the water flow, but instead it seems 
to have a much smaller effect. This paper attempts to 
address the effect of granular properties — the deadzone 
and grain-grain dissipation — upon the cone angle of the 
leaving flow, in order to understand why the granular and 
liquid flows are so broadly similar even though they are 
microscopically distinct. 

Our tool for addressing this problem will be simulation. 
Simulations of the 2D granular jet impact have been per- 
formed using a molecular dynamics approach|12|. The 
grains in these simulations are of flnite rigidity, which 
introduces a fast time scale to the problem (relative to 
the time scale due to relative motion between grains). In 
such simulations, this fast time scale limits the timestep 
used, and therefore the overall computational cost of the 
problem. We instead use a hybrid timestep driven rigid 
body collision method [13] that approximates perfectly 
rigid grains. The perfect rigidity is an approximation 
that we consider supported by the large difference in pres- 



sure scale in the experimental jet (on the order of 100 
kPa, based on glass beads at 10 m/s[T]) and the Young's 
modulus of the grains (on the order of 50 GPa for glass). 
In these simulations we can precisely vary the grain 
properties and the boundary condition in ways that are 
difficult to access experimentally. We expect that the 
cone angles observed in 2D simulations will not quan- 
titatively agree with those observed in 3D experiments. 
However, we are interested in what various granular prop- 
erties do to the cone angle, and we do expect such effects 
to be qualitatively the same in 2D and 3D, as we re- 
produce experimentally measured A — B values with the 
2D simulations. In order to compare to a corresponding 
fluid flow, we can use exact 2D solutions of the Euler 
equations, l4j 



du 
~di 



u-Vu = VP,V 



(1) 



First, we will address the general structure of the flow, 
including the deadzone. Granular material dropped onto 
a surface forms a static pile with a particular character- 
istic angle, the angle of repose, which is determined by 
the properties of the individual grains — their shape and 
the grain-grain friction. We propose that the deadzone 
in granular jet impacts may be a similar feature. Rather 
than being held in place by gravity though, it is held in 
place by the pressure from the incoming granular jet. If 
this is the case, then a necessary element for the forma- 
tion of the deadzone should be whether or not there is 
friction between the grains and the target. If the target 
is frictionless, force chains must come in perpendicular 
to its surface, and a static pile cannot be supported. 

Next, we will address the effects of dissipation in the 
system. There are two sources of dissipation: the coef- 
ficient of restitution (inelasticity), and microscopic fric- 
tion. The restitution coefficient appears to have no effect 
on the cone angle in experiments, and the friction is ex- 
perimentally difficult to control. We use simulations to 
precisely control these parameters and investigate the be- 
havior of the cone angle, and find that there are in fact 
small effects whereby the friction in particular can in- 
fluence the cone angle. We observe inelasticity to have 
almost no effect on the cone angle unless there is abso- 
lutely no friction, and so we do not focus on it too heavily 
here. We propose that friction is playing a mostly per- 
turbative role in this impact geometry, in which case we 
can predict the effect of changing the grain friction by 
computing the total amount of dissipation we expect in 
a homogeneous assemblage of individual grain-grain col- 
lisions. This then acts as an envelope for the amount of 
dissipation experienced by the granular jet, and tells us 
how the cone angle should change. We compare this the- 
oretical prediction to measurements from our simulations 
and find good quantitative agreement. 




FIG. 1: Geometry of the granular jet impact. Grains are 
projected from the left at a target on the right, and form 
a cone with interior angle ^q. In the granular experiments, 
gravity acted in the downwards direction but was negligible 
compared to the velocity scale of the impact. We do not 
simulate the effect of gravity in our simulations. 



CONE ANGLE 

The cone angle may be understood as the combination 
of two factors: transfer of momentum from the target to 
the jet, and dissipation that occurs between the inflow 
and outflow. A derivation for the functional form of the 
cone angle ^o in three dimensions appears inj2j. How- 
ever, care must be taken as the result is different in two 
dimensions. We present a version of the derivation here 
that takes into account an arbitrary dimension d. 

We begin with the flux of momentum parallel to the jet 
(the y direction) per unit time. This is qy ~ QUq, where 
Q is the mass flux Q = pUqctj and Cj is the cross-sectional 
area of the jet and Uq is the inflow speed. The target 
exerts a vertical force on the jet equal to Ft = ApU^at, 
where at is the area of the target and A is a dimensionless 
geometric factor. As such, the y component of the leaving 
velocity is Uy = {qy - Ft)/Q. 

We now use conservation of energy to compute the 
leaving velocity. Some dissipation will occur associated 
with the deformation of the flow around the target. We 
assume that this dissipation is, to first order, propor- 
tional to the area of the target, and so Uf ~ Uo{l~B — ). 
This restricts us to situations in which the target area is 
much smaller than the scale of the overall flow. 

Putting these results together, we have an expression 
for vpo: 



1-A^ 
cos(*o) = Y^7^ 



(2) 



We define x = fr = {jy;)'^~ 
denominator: 



and Taylor expand the 
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FIG. 2: Cone angle measured as a function of Dt/ Dj in 2D 
granular jet impact simulations. Unlike in 3D, the cone angle 
scales as the square root of the ratio of the target to the jet 
diameter. We show data for both stick and slip boundaries, 
at /i = and /x = 1. The dashed black line is the theoretical 
curve for A — B = 0.35. Note that despite the fact that the 
theoretical curve is for the asymptotic limit where x ^ 1, in 
practice the data follow this curve up until x ~ 2. 





FIG. 3: Profiles of the LI norm of the velocity for slip (a) and 
stick (b) boundary condition simulations with /i = 1.0 and 
X = 0.5. In the simulation with a slip boundary condition, 
the velocity linearly approaches zero at the stagnation point, 
whereas in the simulation with the stick boundary condition 
there is a cusp-shaped deadzone within which the velocity 
decays exponentially. The flow is from the top, and the shaded 
square is the target. 
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cos-\l-{A-B)x) 



(3) 



The asymptotic behavior of this equation as x ~^ is 
*o oc {A-B)^ (where x^{%:Y^^)- In 3D, this means 
that the cone angle asymptotes to a linear dependence on 
the ratio of the target diameter to the jet diameter, with 
a slope equal io A — B. In 2D, however, the cone angle 
asymptotes to a square root dependence on the target- 
to-jet ratio, with prefactor A — B. The factor A — B 
describes the combined effect of flow geometry due to 
the target {A) and dissipation (B) and in the limit of 
small target to jet ratio is independent of the size of the 
jet with respect to the target, and so provides a useful 
quantity for characterizing the behavior of the cone angle 
across both 2D and 3D. 

In simulations, we measure this cone angle for x < 
2 by accumulating all grains that have moved past the 
target at least 200 grain radii. We divide the resultant 
grains into those that have x coordinates greater than 
the center line coordinate ('right'), and those that have 
X coordinates less than the center line coordinate ('left'). 
We accumulate the average y component of the velocity 
weighted by grain mass Vy, and the average x component 
of the rate of travel away from the centerline Vx (so grains 



right of the centerline contribute Vx, and grains left of the 
centerline contribute —Vx). In our coordinates, ily < 0. 
We use these velocities to determine the cone angle ^q — 
TT —'ta'n~^ {—Vy / Vx) . When x ^ 2, this technique runs into 
problems as much of the outflow remains above the target 
for the span of our simulation domain. In these cases, we 
must relax the 200 grain radius cutoff, and instead use a 
cutoff that grains should have travelled at least 200 grain 
radii along the direction of the target's surface before we 
consider them. 



We demonstrate the predicted 2D square root scaling 
in Fig. [2] Note that while the scaling only necessarily 
applies in the asymptotic limit as x ~^ 0, it is in prac- 
tice obeyed over almost the entire range of behavior until 
^0 -^ 90°. As such, we may use x — 0.5, which gives us 
a good combination of resolution of flow structures (im- 
proved by a large target) and resolution of the outflow 
cone (best for angles near 45°), and still extract A — B 
from this relation. 
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FIG. 4: Profile of the transverse (a;) velocity along a line 
displaced 0.25-Dt from the centerline. Data are shown for 
two different values of the friction coefficient /x, and for stick 
and slip boundary conditions. In the case of stick boundary 
conditions, the constraint that the velocity be zero at the 
target results in a local maximum in the transverse velocity. 
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FIG. 5: Profile of the longitudinal {y) velocity along the cen- 
terline. Data are shown for two different values of the friction 
coefficient /i, and for stick and slip boundary conditions. In 
the case of stick boundary conditions, the velocity decays ex- 
ponentially towards zero near the target, whereas in the slip 
case it decays linearly. 



DEADZONE 

Experimentally, a granular jet impinging upon a steel 
target produces a deadzone in which the particles are held 
in place by pressure from the incoming jet. In order to 
test whether the deadzone arises from tangential support 
at the target surface, we vary the boundary condition at 
the target in our simulations. We implement two differ- 
ent boundary conditions: 'slip' and 'stick'. In order to 
implement a slip boundary condition, we reverse the per- 
pendicular velocity of grains that impact the target but 
otherwise leave them unchanged. The stick boundary 
condition corresponds physically to a rough target, one 
that is decorated with a layer of grains. We cause any 
grain that impacts the target surface to adhere, and so 
subsequent grains collide with the decorated target, us- 
ing the collision rules for grain on grain collisions. These 
boundary conditions are the two extreme limits of a con- 
tinuum of possibilities. In practice, the experimental tar- 
get has some coefficient of friction for grain-target colli- 
sions, which may be different than that for grain-grain 
collisions. In order to more directly address the matter 
of the presence or absence of the deadzone, we ignore the 
extra degree of freedom here. 

Our simulations consist of a jet one hundred grains 
in diameter, with grains being constantly injected into 
the simulation domain to form the jet, and deleted as 



they leave the simulation domain. When the system is 
fully developed, it contains approximately 140000 grains, 
varying somewhat with cone angle. As this is a system 
comprised of perfect hard spheres, the only dimensionless 
groups are ^ and ^ (where dg is the diameter of a sin- 
gle grain) , and so our results are independent of the jet 
velocity and the grain density. We use a polydisperse set 
of grains uniformly distributed in radius between 0.8 and 
1.2 to prevent crystallization. The grain density is held 
constant, and so the grain mass scales as the square of 
the grain radius. We initialize the jet by placing grains 
randomly in space in an inflow region and then allow- 
ing them to relax their positions to avoid overlap. The 
result is a packing fraction oi (f) — 0.82, close to jam- 
ming. In this paper we will generally discuss results for a 
flxed geometry -^ — 0.5. This is well within the scaling 
regime in which ^q obeys the theoretical prediction of 
Eq. |3] As shown in Fig. [2] we have done simulations at 
several target-jet ratios to determine the quality of this 
scaling assumption, and the data fall precisely upon the 
predicted curve at these values of x- 

For our grain material properties, we choose a coeffi- 
cient of restitution of 0.9, which is consistent with the 
glass beads used in the experimental jet impact. How- 
ever, we also find that our results do not significantly 
depend on our choice of restitution coefficient, as the fre- 
quency of collisions is large enough that almost all avail- 




Ae = tan-i(x 



FIG. 6: Schematic of the set of angles of approach that will 
bring two grains starting at a distance A to collision. 



able energy is dissipated very quickly. We have Coulomb 
friction between the grains, with a coefficient of friction 
/i that we vary. Due to the presence of tangential forces, 
we track both the velocity and angular velocity of grains. 
For the purposes of observing the deadzone, we note that 
whether we obtain a deadzone or not is insensitive to 
our choice of /i, which we have varied from 10~^ to 10. 
This can be seen in Fig. |4] and Fig. [S] which show that 
the general character of the velocity profiles is strongly 
dependent on the boundary condition, but only weakly 
depends on /i. 

We measure the time-averaged velocity profiles near 
the target for both stick and slip boundary conditions. 
We find that in the case of a stick boundary condition, 
the near-target velocity contours have a cusped structure. 
Within that contour, the velocity decays exponentially 
towards zero, with a decay length scale of a few grain 
diameters (Fig. [sb. Fig. [5]). This is consistent with the 
experimentally observed deadzone. However, when we 
switch to a slip boundary condition we observe no dead- 
zone regardless of the grain-grain friction coefficient we 
choose. Instead we observe a stagnation point, where the 
velocity linearly approaches zero (Fig. [sk. Fig. I5|. Sur- 
prisingly, the velocity contours around this stagnation 
point quantitatively agree with the exact solution of the 
2D Euler equation for this geometry, suggesting that even 
despite the high amount of dissipation in the system, the 
flow near the target is dominated by the interaction of 
geometry and the conservation of momentum ^4j. 



FRICTION EFFECTS 

We now turn to the effect of dissipation. We observe 
almost no dependence of the cone angle on inelasticity in 



the system when there is also friction (Fig.^). However, 
the friction appears to have a weak but measurable effect. 
It is this effect that we will try to explain. In principle, 
friction could give rise to a change in deadzone geome- 
try that alters the geometric contribution to the reaction 
force of the target on the jet. However, as we have shown 
in the previous section, this effect is very small (Figs. l4J 
[5]). As such, we will proceed by neglecting the influence 
of friction with deadzone geometry, and observe whether 
this approximation successfully explains the behavior of 
A~B. 

Consider a gas of grains with random velocities. If we 
take two such grains, then the angle between the vector 
between them and their relative velocity will be random. 
We can then restrict ourselves to pairs of grains that will 
collide given their current trajectory. The result is that 
only a narrow band of angles will lead to a collision, but 
that within that band all angles should be equally likely. 
The further apart the grains, the narrower the band of 
angles. This situation is pictured in Fig.lH] The width of 
the angle distribution is a function of the mean distance 
between grains A: 



A9 = tan" 



,2r, 



(4) 



We define the collision parameter b to be the offset 
between the centers of the grain perpendicular to their 
relative velocity at contact, normalized by the sum of the 
grain radii. We can then express the fraction of energy 
lost in a collision as a function of b. We compute the 
distribution of collision parameters in random collisions 
in this granular gas, given an initial distance and uniform 
distribution over colliding values of 6. 

If the two grains approach at an angle 6, then b is: 



where ^ = sin(0) and 



v^d^V^^-\i-(^] e 



(5) 



2r 



2r 



(6) 



Transforming the uniform distribution in 6 to P{b) 
nves us: 



P(6) 



1 



1 



A0\P'i9{b))\ 



(7) 



We show P{b) for various values of A in Fig. [t] 

If we take the limit where A — >■ 2r, then we have a 

situation in which grains are constantly in collision with 

eachother. This corresponds to a dense granular pack. 

While this limit is somewhat questionable, as now there 
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FIG. 7: Collision distribution of in terms of the collision 
parameter b for different average distances between grains. 
When the average distance between grains approaches the 
sum of their radii, the collision distribution diverges with a 
square root singularity near 6=1. Data from a mean-field 
collision simulation using uniformly distributed polydisperse 
grains of radii [0.8, 1.2] is shown for A = 2.1r, and is indistin- 
guishable from the mono-disperse distribution. 



can be a strong ordering in the pack, if the way that 
the pack is driven by exterior forces is random, the basic 
assumption of no correlation between the angle of the 
velocity of a grain and the angle of their contact still 
holds. In this limit, we observe that b — > sin(0), and so 
P{b) asymptotically behaves as: 



P(6)cx 



VT^^ 



(8) 



This result also applies to polydisperse grains in con- 
tact, as the grains do not move between collisions and 
so the entirety of the distribution is determined by the 
random choice of relative velocity, and not the grain ge- 
ometry. This can be seen to hold even away from the 
close-packed limit in Fig. [7J where we present data both 
from analytic calculation of the distribution and from a 
random sampling of 5 x 10^ polydisperse collisions. For 
polydisperse grains, a constant A is not well defined, and 
so instead we hold constant the distance between the 
surfaces of the grains (corresponding to A — 2r in the 
monodisperse case). 

We now look at the energy dissipated in an individual 
collision occurring at collision parameter b, with coeffi- 
cient of restitution e and friction coefficient /i. We as- 
sume that these grains have no relative rotation when 



they come into contact with eachother in order to sim- 
plify the calculation. 

The impulse felt by colliding grains can be separated 
into a normal component In in the direction f and a 
tangential component It in the direction t. If two grains 
approach in the center of mass frame with a momentum 
Pi and no relative rotation, these components are: 



\lN\ = {l + eM-f\ (9) 



|/T|=min(Ml + e)|pl-r|,^) (10) 



The factor I. in \It\ assumes that the grains are homo- 



geneous discs. In general, this factor is , , a/n where 
/ is the moment of inertia of the grain. Observing that 
Pi-r — |pi|"\/l — b"^ and pi-t = \pi\b, we can express these 
impulses in terms of the collision parameter b. Given 
these impulses, we can compute the fraction of energy 
lost in the collision 7 = 1 — -g?^. 

^ -fc-o 

If b/3 < fi{l + e)y/l-b^, then: 



7<(6) = (l-62)(l-6^) + 



562 
IT 



(11) 



otherwise, 



7>(5) = 2(l+e)(l-fe2+^5v/l - b^)-{l+€f{l+p.^){l-b^) 

(12) 
This dissipation function is plotted in Fig. |8] As the 
coefhcient of friction increases, only collisions closer and 
closer to 6 = 1 are affected. However, for a dense pack 
the probability density of such collisions diverges, and 
so these glancing collisions are where the effect of higher 
friction dominates. We estimate the total amount of dis- 
sipation with the integral: 



Dip) = / P{bh{b)db 



(13) 



We expect that this dissipation integral should control 
the dimensionless parameter B in the quantity A — B. 
When /x — >■ 0, we expect A — B to he controlled mostly 
by the geometric contribution (although inelasticity still 
produces some finite dissipation). As such, this asymp- 
totic value is not predicted by the present theory and is 
one fit parameter of this model. Secondly, the dissipa- 
tion integral is an estimate of the fraction of energy dis- 
sipated in a set of uniformly distributed collisions, rather 
than a particular flow. As such, we expect the dissipa- 
tion observed in the jet impact to scale with D but not 
to be exactly equal to D. We therefore have a second 
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FIG. 8: This plot shows the fraction of energy dissipated in 
a colhsion as a function of the colhsion parameter h. 



fit parameter, which is a constant of proportionaUty that 
determines the influence of D upon the quantity A — B. 
As such, we expect that: 
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FIG. 9: a. This plot shows data for A — B from simulations 
at X = 0.5, along with fits to Eq. [14] The value of A — B is 
determined by inverting the equation for the cone angle at a 
fixed value of x- The error bars are determined by measuring 
the cone angle at four different times during the simulation, 
and then computing the standard deviation of the resultant 
A — B values, b. The inset shows the dependence of A — B 
on the coefBcient of restitution for a fixed value oi fi — 0.2. 



A-B^Co-CiD{ii) 



(14) 



To test this, we performed a series of simulations for 
different values of /x, and measured A — B in each case. 
We simulated both shp and stick boundary conditions to 
see if the presence of the deadzone made any difference in 
the scaling oi A—B. The data are shown in Fig. [9^, along 
with fits of the dissipation theory to the data. The values 
of fj, at which A — B switches between its asymptotic 
behaviors are predicted by the theory with no adjustable 
parameters, as the two fitting constants do not influence 
the fj. dependence of the curve, only the vertical scale. 
The first asymptotic regime seems to be fairly insensitive 
to the details of the collision distribution. However, we 
find that if we use a different P{b), e.g. P{b) ex (1 — 6^)°, 
then as a — ^ I the location of the second asymptotic 
regime goes to /i — ?► cxo. 

We obtain agreement within errorbars for both the 
stick and slip boundary conditions. The effect of the 
target boundary condition is a fixed offset between the 
curves — only the fit parameter Co differs between the 
two cases. 



Polydispersity 

We also consider the matter of polydispersity. We gen- 
erally use a polydisperse collection of grains to prevent 
crystallization (which can lead to markedly different be- 
havior in granular systems). Our theory of the cone angle 
seems to be independent of any polydispersity in the sys- 
tem, though it is a mean-field theory of a random pack 
and would not be expected to capture anything depend- 
ing on system-wide ordering. If polydispersity impacts 
the macroscopic behavior of the granular flow, it could 
be an interesting connection between mesoscopic order 
within granular systems (e.g. partial crystallization) and 
bulk behavior. 

Our method for generating the initial jet is to insert 
grains at random in a volume of space well above the 
target. We continuously relax the pack during flight to 
remove the initially generated overlaps. As the jet has a 
free surface, all overlaps are eventually resolved by this 
process (there is no possibility of generating a pack that 
cannot somehow be relaxed). Our grains are generated 
with a uniform distribution of width p, which parameter- 
izes the polydispersity. Normally we use p = 0.4. 

In order to determine whether this polydispersity has 
an influence on A — B, we performed simulations with 
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FIG. 10: This figure shows the correlation function of granular 
packs above the target for polydispersities p = and p = 
0.4. The effect of polydispersity is to decrease the correlation 
decay length. 



a stick boundary condition and /x = 1, but with poly- 
dispersity values 0, 0.1, 0.2, and 0.3. We examine the 
internal structure of these packs above the target (at 
least 400 grain radii up) to identify how the polydisper- 
sity is influencing the microscopic order. We compute a 
two-point correlation function C{r) = - X^ii '^(^ij ^ '')' 
treating each grain as a delta function at its center. The 
result is the correlation function of a solid with disor- 



der (Fig. 10), with several peaks whose amplitudes decay 
exponentially. Measuring the decay constant of the expo- 
nential for each polydispersity we can extract the length 
scale of the ordering. We find that this length scale is at 
most 3.5 grain radii even for the monodisperse pack, and 
decreases to 2.8 for p = 0.3. This suggests that order- 
ing of the initial jet should not play a large role in the 
dynamics of the impact as the length scale does not de- 
pend strongly on the polydispersity. As we increase the 
polydispersity, we decrease the range of the microscopic 
order in our jet. However, even for monodisperse grains 
we never achieve a fully crystalline jet above the target 
due to our method of preparation. 

Because the local crystalline clusters have a different 
length scale than individual grains, we might expect that 
there could be a system size effect associated with the 
polydispersity. We perform simulations for various sys- 
tem sizes for polydispersities oi p — 0, p = 0.2, and 
p = 0.4 (Fig. 11). We find foremost that the monodis- 
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perse case is significantly different than the polydisperse 
cases, with a strong non-monotonic system size effect 



FIG. 11: This figure shows the effects of system size and 
polydispersity for a stick boundary fiow a,t fi — 1. 



with a length scale of Dj/dg — 50. In the monodis- 
perse case, the total range of the system size effect on 
A — B is approximately 0.045, which is comparable in 
scale to the effect of the presence or absence of a dead- 
zone. Furthermore, it is not clear that the monodisperse 
case has reached its asymptote by Dj/dg = 100, and so 
further evolution of the behavior with respect to system 
size may occur. 

For the polydisperse cases. A — B increases with poly- 
dispersity. Aside from the monodisperse case, however, 
A — B is not strongly influenced by polydispersity in the 
large system limit. The scale of the effect we observe 
is approximately 0.004, an order of magnitude smaller 
than the friction and deadzone effects. Even in the poly- 
disperse case, a finite size effect is observed, but with a 
length scale smaller than the monodisperse case (peaking 
at about Dj/dg = 25). The overall impact of system size 
on A — B is also reduced in the polydisperse case, and it 
appears that Dj/dg = 100 is approaching the asymptotic 
limit with respect to system size (in the sense that the 
change oi A — B between a system size of Dj/dg = 50 
and Dj/dg — 100 is about 0.005, and so is about a 2% on 
A — B). In the case of all polydispersities, the A — B de- 
pendence on system size seems to follow the same curve 
for very small systems. At these small values of Dj/dg, 
the outflow is no longer coherent. This suggests that the 
behavior is dominated more by random scattering than 
by bulk flow. As the system size increases, the flow be- 
comes collimated (around Dj/dg = 25), and the curves 
diverge from eachother. Our current model does not pre- 
dict the polydispersity effect at all, but it seems clear 
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that it is a higher-order phenomenon that may not be 
capturable from a simple mean-field approach. 



CONCLUSIONS 

We have shown that the behavior of a granular jet im- 
pinging upon a target is well described by the combina- 
tion of a boundary condition effect at the target to cap- 
ture the geometrical effects of the target reaction force, 
and a microscopic model for the fraction of energy dissi- 
pated in collisions within the jet. The boundary condi- 
tion determines whether or not a deadzone forms, which 
effectively changes the target geometry. A boundary con- 
dition capable of supporting horizontal stresses is needed 
in order to observe a deadzone, whereas a slip bound- 
ary condition does not produce a deadzone even if the 
friction between grains is large. 

The presence or absence of the deadzone is responsible 
for a fixed offset in the value of the constant A— B. On 
the other hand, friction between the grains seems to be 
well-described as a local, homogeneous effect, changing 
only the total dissipation within the jet and not strongly 
influencing the large-scale structure of the flow. While we 
see only a small direct change to the velocity profiles due 
to friction effects (Figs, [i] and |5|, not corresponding to 
any macroscopically altered flow structures, we observe 
a change in cone angle consistent with a microscopically 
homogeneous frictional effect. We have presented a cal- 
culation that estimates the total dissipation within the 
jet as a function of the friction coefficient, and find that 
it predicts the crossover as a function of /i between the 
asymptotic regimes oi A — B corresponding to no friction 
and infinite friction. 

It is not a given that friction should behave in this way. 
In general, factors such as surface tension or viscosity are 
associated with length scales, which give rise to corre- 
sponding extended structures in the fiow (e.g. bound- 
ary layers) along with any homogeneous effects they may 
have. In the theory for the water jet, the infiuence of 
these factors on the cone angle is restricted by going to 
the limit where the length scales associated with surface 
tension and viscosity are much smaller than other length 
scales in the system, and so their effects have asymptot- 
ically saturated |2]. 

In contrast, friction and inelasticity do not create any 



new flow structures and do not introduce any new length 
scales in this flow geometry. As such, the effects are 
well-captured by the homogeneous picture. If we con- 
sider other granular systems, this approach should work 
so long as those systems do not develop boundary layer- 
like structures as a consequence of friction. The success 
of this microscopic description suggests an explanation 
for the correspondence between perfect fluid flow and the 
granular dynamics reported in [13]. The effect of dissi- 
pation is a weak perturbation and does not qualitatively 
affect the flow behavior. 
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